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StttfW  flt©/5  .  ABSTRACT 


Previous  work  on  parameter  plane,  three  dimensional  pa- 
rameter space,  and  singular  lines  in  the  parameter  plane  are 
reviewed. 

A  general  concept  of  n-dimensional  parameter  space  is 
hypothesized  whereby  the  parameter  plane  becomes  a  special 
case  of  the  general  hypothesis.   By  the  same  argument  the 
singular  line  is  shown  to  be  a  special  case  of  the  singular 
hyperplane . 

Existence  criteria  for  singular  lines  are  established, 
and  compensation  methods  for  creating  singular  lines  in  non- 
singular  systems  are  derived  and  used. 
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CHAPTER  I 


INTRODUCTION 


The  analysis  of  linear  control  systems  by  use  of  the 
Laplace  Transforms  of  the  systems1  characteristic  differential 
equations  is  adequately  accomplished  by  root  locus  or  fre- 
quency response  methods  when  all  coefficients  of  the  complex 
variable  s  are  real  constants. 

These  same  methods  are  also  sufficient  to  effect  anal- 
ysis when  the  coefficients  are  functions  of  a  single  real  var- 
iable K.   The  general  characteristic  equation 

n 


f(s)  =  I      a,  sk  =  0  (1-1) 

k=0   K 


may   be    expressed   as 


n  , 

f(s)    =      I       [b    (K)+c    ]sK   =    0    ,  (1-2) 


k=0 

from  which  follows 

I      bk(K)sk 
— =  -1  ,  (1-3) 

I      cks 

k=0 

which  is  then  in  the  proper  form  for  either  root  locus  or  fre- 
quency response  analysis. 

However,  many  modern  control  systems,  although  still  lin- 
early conceived  or  approximated,  cannot  be  handled  by  the 
above  somewhat  "classical"  techniques.   The  complexities 
accompanying  the  advancing  "state  of  the  art"  in  control  system 
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and  compensation  design  often  result  in  characteristic  equation 
coefficients  that  are  functions  of  two  or  more  variables. 

In  1959,  Mitrovic  introduced  a  method  [1]  for  determining 
the  roots  of  the  characteristic  equation  when  the  coefficients 
of  the  two  lowest  ordered  terms  vary  independently. 

A  point  in  the  s-plane, 

s   =   -a   +    joo  (1-4) 

may  be    specified   in  polar   form   as 

s    =    -Coo      +    joj    A-?2.  (1-5) 

n         J    n  ' 

where  co   is  the  radial  distance  from  the  origin  to  the  point 

and  c  is  the  cosine  of  0,  the  angle  between  the  real  axis  and 

00  . 

n 

Since 


and 


£oj   =  cos6  (1-6) 

n 


co  /l-cz  =  sine,  (1-7) 


equation  1-5  may  be  expressed  as 


and 


s  =  a)  (cose  +  jsine)  ;  (1-8) 

n       J 


sk  =  ook(coske+jsinke)  .  (1-9) 

n 


As  previously  defined,  for  left-hand  plane  roots 


therefore 


=  cos"1(-c:)  ;  (1-10) 


coske  =  cosk[cos   (-£)],  (1-11) 


which  can  be  expressed  as  a  Chebyshev  function  of  the  first 
kind,  T,  (-£) ,  and 

sinke  =  sink [cos   (-?)],  (1-12) 

12 


-Re    <■ 


-a 


Fig.    1.1      A  Point   in    the    s-Plane 
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which  can  be  expressed  by  a  Chebyshev  function  of  the  second 
kind,  sin0U,  (-£) .   This  permits  equation  1-9  to  be  written 

sk  =  oj^[Tk(-c)+jsineUk(-c)  ]  (1-13) 

sk  =   cok[(-l)kTk(c)+jsin6(-l)k+1Uk(C)]  (1-14) 

and   equation   1-1   to  be  written 

f(s)    =      I      a,  u)k(-l)kT,  (O+jsinG      I      a.  cok  (-1)  k+1U,  ( ?)    =    0. 
k=0  k=0      K   n  K 

(1-15) 

Since,  for  any  non-zero  value  of  0,  both  the  real  and  imag- 
inary parts  of  equation  1-15  must  independently  equal  zero, 
two  simultaneous  equations  for  a,  are  formed: 

n     .     , 

I      a,/(-l)KT,(c)  =  0  (1-16) 

k=0   K  n      K 

I      a,ajk(-l)k+1U,  (?)  =    0  .  (1-17) 

k=0   K  n        K 

The  following  relationship  exists  between  the  two  kinds  of 

Chebyshev  functions : 


TkU)  =  ^VC)  "  Uk-lU)  '  d-18) 


also 


J       (-l)k+1a.  cokU,   (c)  =    -    I       (-l)ka,  03kU,   (^)  ;  (1-19) 

,  Lk  knks  in  k   n   k 

k=0  k=0 

therefore    for   any  non-zero  value   of    C,    the    simultaneous    equa- 
tions   reduce   to 

X    '-''VnVl'^    =   °  '        (1"20) 

k=0 


k        k. 


I       (-1)    a.u£u.   (C)    =    0    .  (1-21) 

k=0  k   n   k 


14 


U,  (£)  can  be  obtained  for  all  values  of  k  from  the  recur- 
sion formula 

Uk+lU)  =  2?UkU)  "  Uk-lU)'  d-22) 

where,  for  all  values  of  C, 

U_X(U  =  -1  (1-23) 

U0(C)  =   0  (1-24) 

U1(C)  =   1  .  (1-25) 

Equation  1-1  can  be  rewritten  as 

f(s)  =  a  s   +  a   ,s     +  ...  a0s   +  A, s  +  Art  =  0 ,  (1-26) 
n      n-1  2       1      0     ' 

and  equations  1-20  and  1-21  become 

-AQ  +  I     (-l)kak^Uk_l(c)  =  0  (1-27) 

k=2 

n       .    . 

-oj.A.  +  I       (-1)  a,u)*U.  (c)  =  0  .         (1-28) 
n  1    i  o       k  n  k 
k=2 

Thus,  AQ  and  A,  become  the  coordinate  axes  of  the 
Mitrovic  Plane;  and  lines  of  constant  £  and  u  ,  along  with 
constant  real  root  lines  from  the  direct  substitution  of  (-a) 
for  s  in  equation  1-26,  can  be  plotted,  giving  the  solution 
of  the  characteristic  equation  for  any  point  (A0 ,  A,). 

This  method  was  later  extended  by  Elliott,  Thaler,  and 
Heseltine  [2,3]  for  specific  coefficient  pairs,  and  then  gen- 
eralized by  Siljak  [2,4]  for  any  pair  of  variable  coefficients 
in  the  Coefficient  Plane  method,  where  equation  1-1,  when 

written  as 

i  n     v 

f(s)  =  A^1  +  A  sm  +  I      a,sK  =  0,         (1-29) 

1      m     ,  L  _   k 


k=0 
k^l,m 
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results  in  the  set  of  simultaneous  equations  in  (A, ,  A  ) 

^  I'm 


<-i>Vnni-i<s>  +  |-1)VX-iw  +  X  |-1|kaA-i«  =  ° 

k=0 

k^l,m  (1-30) 

1  n  1  m  n  m      ■>  L *  k  n  k  "' 

k=0 

k^l,m  (1-31) 

Although  this  is  an  improvement  of  Mitrovic's  original 
method,  it  is  still  limited  by  the  requirement  that  only  two 
coefficients  of  the  characteristic  equation  contain  variable 
elements.   However,  in  1964  Siljak  [5]  extended  this  process 
into  the  Parameter  Plane,  where  two  variable  parameters  can 
appear  in  any  of  the  coefficients  of  the  characteristic  equa- 
tion.  In  the  linear  case, 

ak  =  bka  +  ck3  +  fk'  (1-32) 

where  a  and  3  are  the  adjustable  parameters.   Hence,  equation 

1-1  becomes 

n     .       n     .     n     . 
f(s)  =  a  I      b.  sK  +  6  I      c,  sK  +  I      f,s   =  0;   (1-33) 
k=0   K       k=0   K     k=0   K 

and  equations  1-20  and  1-21  become 

°  Jo  K,V»UHI(I  +  B  J„  (-1)k°kunUk-l(5' 


+  I      (-l>kfk«>nUk_1<?>  =  0      d-34) 


n 

I 
k=0 


a  J   (-l)kbk^Uk(^  +  e  JQ  ^'VX^ 


k.   k. 


+  I       (-D  fu<U.(5)  =  0  .      (1-35) 
k=0      k  n  K 
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For  notational  ease, 

n 


1  =  I       (-l)kxka)^Uk_1(0  d-36) 


X 

k=0 

X2  =  X    '-"Vn"^1  •  l1"37' 

k=0 

and 

B.oc  +  C13  +  F,  =  0  (1-38) 

B2a  +  C23  +  F2  =  0,  (1-39) 

which  are  solved  by  standard  methods,  giving  constant  £,  ui  , 
and  a  lines  in  the  a-$  Parameter  Plane  for  the  solution  of  the 
characteristic  equation. 

In  1968  Cadena  [6]  extended  the  parameter  plane  concept 
to  three  independent  parameters,  where 

ak  =  bka  +  Ck3  +  dkY  +  fk*  (1-40) 

The  parametric  simultaneous  equations  then  become 

B.a  +  Cie  +  Dxy  +  F1  =  0  (1-41) 

B2a  +  C23  +  D2y  +  F2  =  0,  (1-42) 

which  when  solved  give  the  line  of  intersection  of  two  planes 
in  3-parameter  space.   By  an  ingenious  device  of  rotating  the 
y  coordinate  by  90°  to  form  the  (-3)  ordinate  in  the  a-3  plane, 
Cadena,  by  graphical  projections,  reduces  the  constant  £,,    w  , 
and  a  planes  to  families  of  lines  in  the  0-3  plane  --  each 
line  representing  a  specific  value  of  y.   In  the  same  work,  it 
is  shown  that  for  an  n-parameter  system,  n-roots  may  be  fixed 
by  the  adjustment  of  each  parameter  to  a  specific  value.   In 
this  case  the  remaining  roots  are  then  determined  by  factoring 
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the  quotient  which  results  from  the  division  of  the  character- 
istic equation  by  the  product  of  the  specified  roots. 

What  emerges  from  this  "reduced  characteristic  equation" 
concept  is  that  a  designer  has  at  his  disposal  a  "degree  of 
freedom"  for  each  variable  parameter  in  his  system.   If,  for 
example,  the  system  contains  three  parameters,  three  roots 
may  be  fixed  by  rigid  adjustment  of  all  three  parameters;  if 
only  two  roots  are  specified  (e.g.,  a  complex  conjugate  pair), 
the  reduced  characteristic  equation  then  becomes  a  function  of 
one  parameter  (the  remaining  degree  of  freedom) ,  and  ordinary 
root  locus  methods  determine  the  remaining  roots;  if  only  one 
real  root  is  specified,  the  reduced  characteristic  equation 
is  then  a  two  parameter  case,  leading  to  root  solution  by 
Parameter  Plane  methods. 

Thus  Parameter  Plane  techniques  have  become  a  powerful 
tool  for  both  the  analysis  and  synthesis  of  modern  complex 
control  systems. 
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CHAPTER  II 


THE  SINGULAR  LINE 


In  the  preceding  chapter  it  was  shown  that  constant  5,  u  , 
and  a  lines  could  be  mapped  in  the  a-3  plane  by  the  simulta- 
neous solutions  of  equations  1-38  and  1-39.   Using  Cramer's 
Rule,  points  M  (C,oo  )  can  be  mapped  into  a  point  M(ct,$)  by 


and 


C1F2  -  C2F1 
B1C2  "  B2C1 


(2-1) 


B2F1  -  B1F2 
B1C2  "  B2C1 


(2-2) 


the  point,  M,  being  the  intersection  of  constant  c,    lines  and 

constant  oo   lines.   By  this  method  the  roots  of  characteristic 

n  ■* 

equations  can  be  determined  from  the  output  graphs  of  various 
Parameter  Plane  computer  programs,  using  linear  interpolation, 
if  necessary,  to  find  all  roots  at  any  point  M(a,3). 

An  obvious  limitation  of  this  method  is  that  the  deter- 
minant of  the  coefficient  matrix  be  non-singular,  or  non-zero. 

Although  this  limitation  existed,  it  was  dismissed  as  be- 
ing of  little  significance  until  1967,  when  Bowie  [7],  in 
attempting  to  solve  a  sixth  order  characteristic  equation, 
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found,  that  by  selecting  various  points  M(a,3),  only  certain 
roots  could  be  found  using  Parameter  Plane  methods.   However, 
by  substituting  the  co-ordinates  of  the  point  M(a,3)  into  the 
characteristic  equation  and  factoring,  the  undetermined  roots 
were  found  to  be  complex  conjugate  pairs.   Moreover,  many 
points  M(a,$)  were  found  to  have  the  same  complex  pair  in  com- 
mon; these  common  root  pair  points  formed  a  straight  line  in 

the  a-|3  plane  which  was  the  locus  of  a  constant  ^-constant  co 
r  n 

pair. 

Upon  further  investigation  of  this  phenomenon,  Bowie  dis- 
covered that  lines  of  constant  values  of  L  co  ,  which  he  called 

'   n ' 

singular  lines,  existed  whenever  equations  1-38  and  1-39  were 
linearly  dependent. 

If  the  coefficient  matrix  of  equations  1-38  and  1-39  is 
expanded  in  terms  of  £  and  co  ,  an  infinite  number  of  singular 
lines  can  be  found  by  setting  the  determinant  to  zero  and  solv- 
ing for  either  £  or  oj   in  terms  of  the  other.   For  an  equation 
of  order  n,  there  are  a  maximum  of  2(n-l)  values  of  co  which 
produce  singular  lines  for  each  of  the  infinite  number  of  val- 
ues of  C;  conversely,  (n-1)  values  of  £  will  produce  singular- 
ities for  each  value  of  co  . 

n 

Bowie  illustrated  this  principle  with  two  examples,  which 
will  be  reproduced  at  this  point: 
Example  I 

Consider  the  characteristic  equation 
f(s)  =  0.04s4  +  0.34s3  +  (0.2a+1.12)s2  +  (0 . 5a+$+l . 7) s 

+  23  +  1  =  0.   (2-3) 
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For   ease    of    coefficient    identification,    display    the    terms    as 


- 

- 

0.2 

0.5 

- 

- 

- 

- 

1 

2 

0.04 

0.34 

1.12 

1.7 

1 

Then ,    if    C   =    0.5, 


B,    =    0.2oj 
1  n 


C1   =    -2 


B~    =    -0.5oo      +    0.2co  C~    =    -co       ; 

2  n  n  2  n 


(2-4) 


and 


B,C0    -    B0C,    =    -0.2co3    +    0.4co2    -    go      =    0 . 
1    z  A   l  n  n  n 


(2-5) 


Reducing , 


to      -    2co      +5    =    0; 
n  n 


(2-6) 


and 


oj      =    1    ±    j2    . 
n  J 


(2-7) 


Since  no  real  values  of  co   exist,  no  singular  line  will  exist 
for  this  characteristic  equation  for  the  value  £  =  0.5. 
Example  II 

f(s)  =  s6  +  80s5  +  (20a+1600)s4  +  840as3  +  (1600a+400B) s2 

+  16003s  +  16003  =  0.   (2-8) 


21 


Displaying  terms, 


- 

- 

20 

840 

1600 

- 

- 

- 

- 

- 

- 

400 

1600 

1600 

1 

80 

1600 

- 

- 

- 

- 

For    C   =    0.5, 


B,    =    1600oo2    -    840oo3 
1  n  n 


B„    =    1600co2    -    20oo4 
2  n  n 


C,    =    -1600    +    400oo 
1  n 


C0    =    -I6OO00      +    400oj 
2  n  n 


(2-9) 


and 


Solving, 


oj4    -    42oo3   +    164oo2    -    320oj      +    320    =    0. 
n  n  n  n 


00      =    2.1115,    37.889,    and   a   complex  pair. 


(2-10) 


(2-11) 


Therefore,    for    c   =    0.5,    there   exist   two   values   of   oj    — 2.1115 

'  n 

and  37.889,  which  will  produce  singular  lines  for  equation  2-8 

Once  the  conditions  for  singularity  have  been  met,  the 
singular  line  may  be  represented  by  the  equation 


3  =  - 


Bl     Fl  _    B2     F2 


(2-12) 


"1     ^1      "2      2 

Cadena  [6]  investigated  the  singular  line  when  the  char- 
acteristic equation  coefficients  are  functions  of  three  var- 
iables . 

Since  the  solution  of  equations  1-41  and  1-42  is  a  line 
in  three  dimensional  space  for  a  C-00  pair,  the  definition 
singular  line  may  be  loosely  applied.   Where  that  line  inter- 
sects the  a-3-co-ordinate  plane  is  the  normal  two  dimensional 
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solution  (or  singular  point)  of  equations  1-38  and  1-39,  since 
Y  =  0.   This  actually  reduces  to  his  "degree  of  freedom"  con- 
cept, whereby  two  roots  may  be  specified  by  specifying  the 
proper  values  of  two  of  the  variable  parameters. 

A  special  case  of  the  three  dimensional  singular  line  is 
when  it  happens  to  lie  in  one  of  the  co-ordinate  planes  or  in 
a  plane  parallel  to  one  of  the  co-ordinate  planes.   When  this 
occurs,  the  singular  line  of  equation  2-12  exists  for  the 
specific  value  of  y • 

An  illustrative  example  of  the  foregoing  was  given  by 
Cadena: 
Example  III 
Given 
f(s)  =  s4  +  (10a+10)s3  +  (40a+5B+10y+30)s2  +  ( 80a+103+252) s 

+  70a  +  103  +  25y  +  340  =  0,  (2-13) 
it  is  desired  that  a  pair  of  complex  conjugate  roots  be  lo- 
cated at  c  =  0.5,  a)   =2.   Division  of  equation  2-13  by  the 

2 
product  of  the  desired  roots,  s   +  2s  +  4 ,  yields  a  reduced 

characteristic  equation  in  a ,  3,  y  and  a  remainder  which  must 
equal  zero  for  the  division  to  be  exact. 

(-20y+200)s  +  (-10a-103-15y+300)  =  0        (2-14) 
or 

R1S  +  R0  =  °  '  (2-15) 

Obviously,  R,  ,  and  RQ  must  simultaneously  be  zero,  so 

y  -  10  =  0  (2-16) 

and 

10a  +  103  +  15y  -  300  =  0.  (2-17) 
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This  yields  the  singular  line 

a  +  3  -  15  =  0  (2-18) 

in  the  plane  y   =   10. 

Perhaps  a  more  illuminating  manner  in  which  to  view  this 
problem  is  to  display  the  coefficients  in  tabular  form  as  was 
done  for  the  two-dimensional  examples.   Thus 

k    4    3    2     1     0 


- 

10 

40 

80 

70 

- 

- 

5 

10 

10 

- 

- 

10 

- 

25 

1 

10 

30 

252 

340 

Then,  for  C  =  0 . 5 ,  oo   =2, 

i  ^  n 


Bl  =  10    Cl  =  10    Dl  =  15    Fl  =  "300 
B2  =   0    C2  =   0    D2  =  40    F2  =  -400 


(2-19) 


and 

10a  +  103  +  15y  -  300  =  0  (2-20) 

40y  -  400  =  0.  (2-21) 
This  simultaneous  system  is  the  singular  line  formed  by  the 

intersection  of  the  two  planes  described  by  equations  2-20  and 
2-21.   The  singular  line  is 

a  +  3  -  15  =  0  (2-22) 
in  the  plane 

Y  =  10.  (2-23) 
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Now,  if  dL  =  10, 


D2  =  20;  (2-24) 


and 

10a  +  106  +  15y  -  300  =  0  (2-25) 

20y  -  400  =  0.  (2-26) 

The  singular  line  is  then 

a  +  6  =  0  (2-27) 

in  the  plane 

Y  =  20.  (2-28) 

If  now  d,  =  d3  =  d,  =  10 , 

D   =95     D2  =  180;  (2-29) 

10a  +  106  +  95y  -  300  =  0  (2-30) 

180y  -  400  =  0;  (2-31) 

and  the  singular  line  becomes 

a  +  6  -  ^   =  0  (2-32) 


in  the  plane 


Y  =  ^S-  .  (2-33) 


9  ' 
By  now  it  is  apparent  that  any  variation  of  coefficients  d, 

will  result  in  different  singular  lines  in  different  planes 

parallel  to  the  a-6  plane. 

What  then  happens  if  coefficients  b,  and/or  c,  are  varied? 

Let  c3  =  c4  =  10.   Then 

C1   =  -70     C2  =  160  (2-34) 

and 

10a  -  706  +  15y  -  300  =  0  (2-35) 

-1606  +  40y  -  400  =  0.  (2-36) 
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The  singular  line  for  £  =  0 .  5 ,  co   =2  now  exists  in  a ,  3;  Y~ 
space  at  the  intersection  of  the  two  planes  described  by  equa- 
tions 2-35  and  2-36.   Of  what  use  is  this  information?   If 
equation  2-36  is  solved  for  y   with  the  result  then  substituted 
into  equation  2-35,  the  singular  line  is  then  "mapped"  into 
the  a-3  plane  as 

a  -  3  -  15  =  0  (2-37) 

via  the  "transform" 

y  =  43  +  10.  (2-38) 

What  does  this  mean?   For  any  consistent  three  dimensional  set 
of  equations  1-41  and  1-42,  a  two  dimensional  singular  line 
can  be  realized  by  placing  an  appropriate  constraint  upon  the 
third  parameter. 

For  further  examination,  let  c~  =  0.   Then 

C,  =  -10     C2  =  -20;  (2-39) 

10a  -  103  +  15y  -  300  =  0  (2-40) 

-  203  +  40y  -  400  =  0;  (2-41) 

4a  -  3  -  60  =  0  (2-42) 

Y  =  y3  +  10.  (2-43) 


and 


if 


Let  b.  =  10.   Then 


B2  =  160;  (2-44) 


10a  +  103  +  15y  -  300  =  0  (2-45) 

160a  +  40y  -  400  =  0;  (2-46) 
and 

50a  -  103  +  150  =  0  (2-47) 
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if 

Y  =  -4a  +  10.  (2-48) 

What  has  been  demonstrated  in  the  preceding  examples  is 
the  existence  of  a  powerful  tool  for  the  synthesis  of  self- 
adaptive  systems  or,  for  that  matter,  any  system  for  which  a 
characteristic  response  is  desired.   How  this  tool  is  used 
will  be  the  subject  of  further  discussion  later  in  this  thesis 
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CHAPTER  III 

N-DIMENSIONAL  PARAMETER  SPACE 

Consider  a  set  of  n  consistent,  independent  simultaneous 
equations  where  n  =  N: 

B-a  +  C-,6  +  D   y   +  E,  6  +  .  .  .  +  N..V  +  P.,  =  0 

B2a  +  C23  +  D2y  +  E26  +  ...  +  N2v  +  P2  =  0 

(3-1) 

Ba  +  C6  +  Dv  +  E(S+...+Nv  +  P   =0. 
n     n     n'     n  n     n 

The  solution  of  the  system  of  equations  3-1  is  the  set  of 
values  obtained  for  the  variables  a,  3,  y,    6,  ...,  v,  which 
defines  a  point,  which  is  a  zero  (n-n)  dimensional  entity,  in 
n-dimensional  space. 

If  now  only  (n-1)  independent  equations  in  n  =  N  var- 
iables can  be  obtained,  the  set  is  reduced  to 

B,  a  +  C,  3  +  D,  y  +  E,6  +  .  .  .  +  N;,v  +  P   =  0 

B2a  +  C26  +  D2y  +  E26  +  ...  +  N2v  +  P2  =  0 

(3-2) 

B   ,a  +  C   n6  +  D   ,y  +  E   -6  +  ...+N   nv  +  P   n  =  0. 
n-1     n-1     n-1'     n-1  n-1     n-1 

The  solution  to  set  3-2  is  a  one  [n-(n-l)]  dimensional  line 

defined  by  an  arbitrary  pair  of  the  N  variables. 

Similarly,  for  (n-2)  independent  equations,  the  solution 

becomes  a  two  [n-(n-2)J  dimensional  plane  defined  by  an 

arbitrary  trio  of  the  N  variables. 
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By  induction  it  is  clear  that  if  only  two  independent 
equations  describe  a  system  of  N  variables,  there  results  equa- 
tions 

B.,a  +  C,B  +  DjY  +  Ey6  +  ...  +  NjV  +  P.=  0     (3-3) 

B2a  +  C23  +  D2y  +  E26  +  ...  +  N2v  +  P2  =  0     (3-4) 

whose  solution  is  an  (n-2)  dimensional  hyperplane  defined  by 
an  arbitrary  (n-1)  set  of  the  N  variables. 

If  the  coefficients  B,  C,  D,  E,  .  ..,  N,  P  are  now  func- 
tions of  some  constant  pair  (C,co  ),  the  solution  to  equations 
3-3  and  3-4  is  an  (n-2)  dimensional  hyperplane  of  constant  £ , 

03  .   If.  when  all  but  two  arbitrary  variables  are  set  to  zero 

n  2 

and  equations  3-3  and  3-4  remain  consistent  and  independent, 
the  point  of  intersection  of  the  (a-2)  dimensional  constant  c; , 
co   hyperplane  with  the  plane  of  the  two  remaining  co-ordinates 
is  determined.   If  equations  3-3  and  3-4  become  dependent,  the 
intersection  is  a  line;  if  consistent,  no  intersection  exists. 

For  non-linear  combinations  of  the  N  variables,  the  (n-2) 
dimensional  hyperplane  then  becomes  hyperbolic,  parabolic, 
cubic,  quartic,  etc.,  in  nature,  which,  in  the  independent  and 
consistent  case  described  in  the  preceding  paraqraph,  gives 
rise  to  the  possibility  of  multiple  intersections  with  the  two- 
dimensional  co-ordinate  planes. 

Although  extremely  little  of  practical  engineering  value 
can  be  extracted  from  the  preceding  discussion,  the  concepts 
expressed  are  valid  and  are  useful  for  an  appreciation  of  the 
two  parameter  problem.   Since  mental  imagery  extends  at  best 
to  a  relatively  poor  perception  of  three  dimensions,  the  root 
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finding  problem  is  limited  at  the  present  state  of  the  art  to 
the  three  dimensional  case  demonstrated  by  Cadena  [6]. 

However,  one  limited  use  is  the  determination  of  necessary 
values  of  two  parameters,  say  a  and  3/  given  arbitrary  values 
of  the  remaining  n  parameters,  for  obtaining  a  desired  £ ,  co 
pair. 

Given  a  system  described  by  the  characteristic  equation 

n     . 
f(s)  =  I      a.  s   =  0  (3-5) 

k=0   K 

where 

ak  =  bka  +  Ck^  +  dkY  +  ek5  +  **•  +  pk'       (3-6) 
one  arrives  at  equations  3-3  and  3-4,  where  all  capital  let- 
ter coefficients  are  functions  of  the  desired  t.    co   pair.   If 

'   n  " 

variables  y,    6,  ...,  v  can  be  fixed  or  measured,  the  problem 
then  becomes  two  dimensional  in  a  and  3; leading  to  solution 
by  the  methods  outlined  in  chapter  one.   This  is  essentially 
what  Cadena  achieves  by  graphical  means  by  folding  the  y~co- 
ordinate  into  the  a-3  plane. 

Figure  3.1  is  illustrative  of  the  three  dimensional  linear 
case.   For  the  system 

Bxa  +  Cx3  +  Dxy  +  F1   =  0  (3-7) 

B2a  +  C23  +  D2y  +  F2  =  0,  (3-8) 

the  constant  l,  ,  co  [{n=3}-2]  hyperplane  is  a  line  in  three 
dimensional  space. 

For  the  case  of  non-linear  combinations  of  variables, 
the  C/  w   lines  become  curves  in  space  with  the  possibility 
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,w   with 
n 

with  plant 


> 


Fig.  3.1   Intersection  of  Two  Planes  in  Three  Parameter 
Space 

31 


of  multiple  intersections  with  constant  y-planes  as  illustrated 
in  figure  3.2  for  the  system 

B^  +  C^  +  D-y  +  H1  =  0  (3-9) 

B2a  +  C26  +  D2y  +  E2a2  +  F^2    +  G2y2  +  H2  =  0.   (3-10) 

For  the  condition  when  equations  3-3  and  3-4  are  depend- 
ent, the  solution  is  an  (n-1)  dimensional  hyperplane  described 
by  either  of  the  two  "singular"  equations.   The  same  line  of 
reasoning  applies  as  in  the  "non-singular"  case,  resulting  in 
the  singular  plane  for  the  three  parameter  problem  and  the 
singular  line  for  the  two  parameter  problem,  which  is  illus- 
trated in  figure  3.3. 

Although  the  singular  line  is  the  most  useful  of  all  the 
"singular"  hyperplanes  because  of  the  availability  of  graphical 
means  for  determining  the  remaining  roots  of  a  characteristic 
equation,  singular  theory  can  be  readily  applied  to  any  multi- 
variate self-adaptive  system  or  any  other  system  where  a  con- 
stant complex  root  pair  is  desired.   For  any  system  with  n 
parameters,  variation  in  any  or  all  of  (n-1)  parameters  may  be 
compensated  by  adjustment  of  the  nth  parameter. 
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Curve  of  constant 
C,u)  with  points 

of  intersection 
with  plane  y=K 


*2a+C23+D2Y+E2a2+F232+G2Y2+H2=0 


B1a+C;,6+D1Y+H1=0 


Fig.  3.2  Non-linear  Intersection  in  Three  Parameter  Space 
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Singular  line 
B1a+C16+F1=0 


Singular  line 
C16+D1y+F1=0 


a 


Singular  plane 
B1a+C16+D1Y+F1=0 


Singular  line 
Bla+D1Y+Fl=0 


A 


Y 


Fig.  3.3  Singular  Plane 
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CHAPTER  IV 

INTRODUCTION  TO  COMPENSATION 

Since  the  singular  line  can  be  an  extremely  useful  system 
characteristic,  the  ability  to  design  compensation  to  produce 
a  singular  line  for  a  specific  complex  conjuqate  root  pair  be- 
comes very  desirable.   The  simplest  mathematical  scheme  which 
comes  to  mind  is  that,  for  a  system 

B,a  +  C13  +  F1  =  0  (4-1) 

B2a  +  C2B  +  F2  =  0 ,  (4-2) 

a  singular  line  is  present  if  B,  =  B?  ,  C,  =  C2,  and  F,  =  F2 . 
Using  the  general  coefficient  X,  if 

Xx  =  X2,  (4-3) 


then 


and 


Xl  "  X2  =  0/  (4"4) 


n      .    . 

I       (-1)  \v\.i(S}-Uk(C)]  =  0,  (4-5) 

k=0 

which,  when  expanded,  becomes 

2  2   3 

-  xQ  +  conx1  +  (l-2c)conx2  +  (-1-2^+4?  )  wnx3 

+  (-l+4c+4c2-8c3)o)4x4  +  (l+4c-12c2-8<;3+16?4)^x5 

+  (l-6^-12c2+32c3+16?4-32c5)oo6x6 

*  (-l-6?+24c2+32c3-80C4-32c5+64,6)^x7 

+  (-l+8c+24<;2-80c3-80c4  +  192C5+64c6-12  8c7)oo8x8 

+  (1+8c-40?2-80C3+240c4+192?5-44  8?6-12  8c7+2  56c8)oj9xq 

n    y 

+    (l-10c-4  0C2+16  0<;3+240S4-6  72£5-44  8c6  +  1024c7+2  56c8 

-512c9)oo10xin    +    ...    =    0    .  (4-6) 

n      10 
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If  the  coefficients  of  X>  for  each  term  k  are  displayed 
as  in  Table  4-1,  the  result,  neglecting  the  k  =  0  row,  which 
is  (-1)  for  m  =  0,  is  a  lower  left  diagonal  matrix  with  each 
coefficient  repeating  once  in  the  same  column.  The  repeating 
coefficient  pairs  alternate  in  sign  as  each  column  is  formed. 
The  values  of  the  main  diagonal-pair  are  (-2)  ,  and  the  suc- 
ceeding diagonal-pairs  are  formed  by  entering  Table  4-2  with 
arguments  1  and  m  to  obtain  the  necessary  multipliers  of  the 
magnitude  of  the  elements  in  the  main  diagonal-pair. 

Table  4-2  is  easily  formed  by  first  making  elements  ln , 

m .  and  1^ ,  mn  equal  to  one ;  then 

1.,  m.  =  1.,  m.  n  +  1.  ,,  m.  .  (4-7) 

l    j  l    j-1  l-l    3 

Therefore,  by  establishing  the  preceding  conditions  for 
all  coefficients  X  of  equations  4-1  and  4-2,  compensation  to 
produce  a  singular  line  can  be  effected. 

Because  plant  gain  is  usually  the  available  control  to 

adjust,  it  is  convenient  to  equate  all  terms  x„  in  the  form 

xn  =  co  x,  +  (1-2O0J  x„  +  ...  ,  (4-8) 

0     n  1  n  2 

since  gain  for  a  type-one  unity-feedback  plant  with  no  zeros, 

which  can  always  be  realized  mathematically  by  block  diagram 

algebra,  is  the  sum  of  all  the  x~  terms. 

This  compensation  can  be  easily  achieved  with  a  cascaded 

component  whose  transfer  function  numerator  contains  the  term 

bna  +  cn3  +  fn  and  whose  denominator  cancels  the  numerator  of 

the  uncompensated  plant. 
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Table  of  Multipliers 
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Example  IV 

Example  I  in  chapter  II  presents  a  characteristic  equation 

0.04s4  +  0.34s3  +  (0.2a+1.12)s2  +  (0 . 5a+3+l . 7) s  +  23  +  1  =  0 

(4-9) 
which  has  been  shown  to  contain  no  singular  lines.   Suppose, 
for  example,  it  is  desired  to  have  a  singular  line  for  the 
roots  £  =  0.5,  go   =  1.   If  equation  4-9  represents  the  system 
in  figure  4.1a,  addition  of  a  cascade  compensator  with  trans- 
fer function 

G  (S)  =  ba+cS+f  (4-10) 

^c{S>  23+1  {q        ' 

results  in  a  new  characteristic  equation 

0.04s4  +  0.34s3  +  (0.2a+1.12) s2  +  (0 . 5a+3+l . 7) s 

+  ba  +c3  +  f  =  0,      (4-H) 

where,  if  the  system  is  made  singular  by  definition, 

B1   =  -b  +  0.2  =  B2  =  -0.3  (4-12) 

C±   =    -c  =  C2  =  -1  (4-13) 

F,  =  -f  +  0.78  =  F2  =  -0.62;  (4-14) 

and  the  system  then  has  a  singular  line  solution  for  the  roots 

C  =  0 . 5 ,  oo   =1  when  b  =  0.5,  c  =  1,  and  f  =  1.4.   The  com- 
'   n  '       ' 

pensated  system  is  then  as  shown  in  figure  4.1b;  and  the  sin- 
gular line  for  the  desired  roots  is 

0.3a  +  3  +  0.62  =  0.  (4-15) 

Although,  as  stated  in  chapter  two,  a  system  may  possess 
an  infinite  number  of  singular  lines,  one  for  a  desired  root 
pair  may  still  not  exist  without  compensation. 
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Fig.  4.1   Example  IV-4th  Order  Plant 
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Example  V 

Figure  4.2a  shows  a  system  which  was  contrived  to  have 
singular  lines.   The  characteristic  equation  is 


s3  +  as2  +  (K1+8K1K23)s  +  20KXK23  =  0, 


(4-16) 


where 


B,    =   oo 
1  n 


B      =    2£u> 
2  n 


C,    =    -20K,Ko         C0    =    -8KnK0(ja 
1  12  2  1    2    n 


F,    =    -2^00 
1  n 


3  2    3 

F~    =    oo      -    4c    oo      -   K-.0J    . 
2  n  n  In 


(4-17) 
(4-18) 
(4-19) 


By  the  definition  for  a  singular  line, 


Bi  ci 


B2   C2 


-F    C 
*1   ^1 


-F    C 
*2   U2 


Bl   "Fl 


B2   -F2 


=  0 


(4-20) 


Substituting , 
2 


00 


n 


-20K1K2 


2?w  -8KnK„ 

n  12 


=    40K,KoCw2    -    8K-.K„oo2    =    0  (4-21) 

1    2.      n  1    2    n 


2£u) 


n 


"20K1K2 


(4c2oo3-w3+Knoo    )         -8KnK0 
n      n      1   n  12 


=    80K,K0£2oo3    -    20KnK„oo 
12^      n 


3 

l"2~n 


+    20K2K„oj      -    16K,K~Coj3         n 
12n  1    2      n   =    0 


(4-22) 


CO 


n 


2<;oo 


n 


2?oo2       (4  6:2oo3-oo3+Knoo    ) 
'    n  n      n      1   n 


=    4c2oo5    -    oo5    +   Knoo3    -    4c2oo3    =    0. 
s      n  n  In  n 


Solving,  equation  4-21  gives 

C  =  0.2, 


(4-23) 


(4-24) 


and  equations  4-22  and  4-23  give 


oo   —  K,  . 
n    1 


(4-25) 
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The  singular  line  for  roots  t,    -    0.2,  oo   =  K-  is  described  by 
both 

a  -  20K23  -  0.4/K^  =  °  (4-26) 

from  the  B, ,  C,  ,  F,  coefficients  and  by 

0.4a  -  8K2/K^3  -  0.16^  =  0  (4-27) 

from   the   B2 /    C2 ,    F2    coefficients.      Dividing   equation    4-27    by 
0.4, 

a    -    20K2/K7g    -    0.4/kT  =    0.  (4-28) 

Since   equations    4-26    and   4-28    must  be    identical,    it   is   ob- 
vious   that 

/K~  =   k.    =   go2   =   u>      =1.  (4-29) 

linn 

This  is  fortunate  if  the  roots  t   =  0 . 2 .  co  =1  are  de- 

'   n 

sired,  but  suppose  a  singular  line  is  wanted  for  the  roots 

5  =  0.5,  oo   =1. 
'   n 

By  block  diagram  algebra  the  system  can  be  reduced  to 

that  shown  in  figure  4.2b.   Selecting  a  cascade  compensator 

which  cancels  the  plant  numerator 

8K,K9gs+ba+cg+f 
Gc(s)  =  SKJ^s+i.S)  (4-30) 

where  the  s   numerator  coefficient  must  contain  $  in  order 
that  C?  be  non-zero  and  where  8K,K2  is  inserted  with  the  idea 
of  hopefully  preserving  a  resemblance  of  the  original  gain,  a 
new  characteristic  equation 

s3  +  as2  +  (K1+8K1K26)s  +  ba  +  eg  +  f  =  0     (4-31) 

is  formed.   For  singularity 
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B,  =  -b+1  =  B2  =  1;         b  =  0  (4-32) 

C1  =  -c    =  C2  =  -8K1K2;   c  =  SK^  (4-33) 

F1  =  -f-1  =  F2  =  -K1;      f  =  K1-l  .         (4-34) 

If  now  K1  =  K2  =  1,  the  compensator  becomes  the  readily  real- 
izable 


_s+l 
c  "~     s- 


G~(s)  =  ^33  '  (4"35) 


with  the  compensated  system  in  figure  4.2c  having  the  singular 
line 

a  -  83  -  1  =  0.  (4-36) 

Figure  4.3  shows  the  parameter  plane  plot  of  the  charac- 
teristic equation 

s3  +  as2  +  (80+l)s  +  86  =  0  (4-37) 

with  the  singular  line  of  equation  4-36 .   The  intersection  of 
the  singular  line  with  the  real  root  lines  then  determines 
the  complete  solution  of  equation  4-37  for  selected  singular 
line  values  of  a  and  3.   It  is  interesting  to  note  that  the 
plots  of  the  constant  £  and  constant  go   lines  form  saddles  as 
they  approach  the  singular  values.   This  is  due  to  computer 
truncation  as  the  program  solution  approaches  the  indeterminate 
condition . 
Example  VI 

For  the  plant  shown  in  figure  4.4a  with  characteristic 
equation 

s3  +  (a+2)s2  +  2as  +  KK-  =  0,  (4-38) 

cascade  compensation  is  desired  to  produce  a  singular  line 

for  roots  £  =  0 . 5 ,  go   =1.   The  first  step  is  to  cancel  the 

n  r 

43 


Scale:  a=5/in. ,&=l/in. 
Fig.  4.3  Example  V-Charactenstic  Equation 
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plant  numerator  while  concurrently  introducing  the  second  pa- 
rameter in  at  least  two  different  characteristic  equation  co- 
efficients so  that  neither  C,  nor  C~  are  zero.   This  can  be 
accomplished  in  at  least  two  ways : 

r    fe\    -    3s+ba+cg+f  ..    ,  v 

GC(S)  -  KKl(s+K2)  (4"39) 


or 


where  the  a3  term  in  the  numerator  is  needed  because  of  the 
a3  term  produced  by  the  denominator.  The  first  method  pro- 
duces the  compensated  characteristic  equation 

s4  +  (a+K„+2)s3  +  [ (K0+2)a+2K0]s2  +  (2K0a+3)s  +  ba  +cB  +  f  =  0, 


with 


(4-41) 


B1  =  -b  +  1  +  K2  =  B2  =  2  -  K2  (4-42) 

C1  =  -c  =  C2  =  -1  (4-43) 

F1  =  -f  +  K2  -  2  =  F2  =  2K2  -  1  ,         (4-44) 


resulting  in  the  singular  line 


(2-K2)a  +  3  +  1  -  2K2  =  0.  (4-45) 

The  compensator  is  then 

3{s+l+[(2K9-l)a-K9-l]/$} 
<VS>  =  KKl(s+K2)   •         (4"46) 

which  is  shown  in  figure  4.4b. 

With  the  second  method  the  characteristic  equation  becomes 

s4  +  (a+3+2)s3  +  (2a+23+a3)s2  +  (2a3+l)s  +  ba  +  eg 

+  da3  +  f  =  0,    (4-47) 
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where 


making 


B1   =  -b+1  =  B2  =   2  (4-48) 

C,  =  -c+1  =  C2  =   2  (4-49) 

D1  =  -d+1  =  D2  =  -1  (4-50) 

F,  =  -f-2  =  F2  =  -2,  (4-51) 


_  ,    v    s-a-3+2a3  #,■  Co\ 

Gc(s)  =  KKl(s+e)     '  (4"52) 


which  is  shown  in  figure  4.4c.   The  singular  line  is 

2a  +  20  -  a3  -  2  =  0,  (4-53) 

which  may  be  written 

3  =  -£*•  +  2  (4-54) 

a-2 

and  is  plotted  in  figure  4.5.   Since 

lim  3=2,  (4-55) 

a->oo 

3  could  be  set  to  2  to  give  singular  line  performance  for 
small  variations  about  a  large  value  of  a. 

The  preceding  has  demonstrated  a  rather  simple  mathemat- 
ical approach  to  singular  line  compensation,  but  it  is  suffi- 
cient both  for  a  practical  solution  to  the  singular  line 
problem  and  for  an  appreciation  of  a  more  elegant  method  which 
is  the  subject  of  the  next  chapter. 
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CHAPTER  V 

COMPENSATING  TECHNIQUES 

Consider  the  uncompensated  two  parameter  control  system 
described  by 

B^a  +  C.^  +  F,  =  0  (5-1) 

B2a  +  C23  +  F2  =  0.  (5-2) 

If  this  system  is  now  compensated  by  some  function  of  s  and 

Y,  where  y  is  a  function  of  a  and  $,  equations  5-1  and  5-2  be- 
come 

Bxa  +  C13  +  D1y(a/3)  +  F,  =  0  (5-3) 

B2a  +  C2B  +  D2y(a,B)  +  F2  =  0  (5-4) 

for  the  compensated  system.   After  equation  5-4  is  multiplied 
by  D^/D2,  the  set  becomes 

Bxa  +  C13  +  D1Y(a/6)  +  F±   =  0  (5-5) 

D  B       DC  D  F 

-jj-^-  a  +  -~^-   3  +  DlY(a,3)  +  -i-=-  =  0.       (5-6) 

Adding  equations  5-5  and  5-6  gives 

D  B  DC  D  F 

(B1  +  l5p0  +  (ci  +  "5^*  +  2V(a'3)  +  Fl  t  -kf  =   °' 

(5-7) 
which,  when  solved  for  y(«»3)»  yields 

1   Bl    B2       Cl    C2      Fl    F2 
Y(o,B)  =  .  l[(_i+_i)a+  (gi+^B  +  5i+  of]'*5-8' 

the  value  of  y  f°r  the  compensator.   Also, 


49 


n       /       f»          ,"D2B1-D1B2,             #-D2Dl"DlC2%0         °2F1         D1F2 
DlY(a,3)     =    ( 2D^ >a    +    < 2D^ >3    "    2DJ-   "    2DJ"    ' 

(5-9) 

which,  when  substituted  into  equations  5-5  and  5-6,  gives 

D2B     D  B  D2C     DC  BjF, 

<B1  -  25^  -  2^>  °  +  (C1  "  25^  "  3^  B  +  Fl  "  ^ 

-   25^  "  °       <5"10> 

,D1B2    D2B1    D1B2,    ,  ,D1C2    D2C1    D1C2,  .  ,  D1F2 
(  D2   "  2D2   "  2D2   a      D2  "  2D2   "  2D2  '      D2 

DP     D  F 

Reducing,  equations  5-10  and  5-11  become 

(D2B1-D1B2)a  +  (D2C1-D1C2)3  +  ^2F1    "  °1F2  =  °   (5_12) 
(D1B2-D2B1)a  +  (D1D2-D2C1)6  +  D^  -  D^  =  0,  (5-13) 

which,  differing  by  a  constant  factor  of  (-1) ,  are  obviously 
singular  and  are  indeed  the  equation  of  a  singular  line  for 
the  compensated  system. 

Several  examples  utilizing  equation  5-8  to  design  a  com- 
pensator to  produce  a  singular  line  described  by  equations 
5-12  and  5-13  will  be  given  for  a  variety  of  systems. 
Example  VII 

Type  Q-3rd  Order-Feedback  Compensation-One  Parameter  in  Plant 
Denominator 

For  the  system  shown  in  figure  5.1a  it  is  desired  to  obtain 
feedback  compensation  to  produce  a  singular  line  for  the  roots 
C  =  0 . 5 ,  oj   =1.   The  first  step  is  obviously  to  introduce  the 
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Fig.  5.1  Example  VII-Type  0,  3rd  Order  Plant 
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second  parameter.  A  slight  amount  of  caution  must  be  exercisec 
to  insure  that  the  second  parameter  occurs  in  coefficients 
other  than  those  which  will  contain  the  compensating  function 
y(a,3)  so  that  the  3  coefficient  in  the  singular  line  equatio 
5-12  will  be  non-zero.  This  is  shown  in  figures  5.1b  and  5.1c 
which  features  the  introduction  of  the  compensating  function. 
The  resultant  characteristic  equation  is  then 

s3  +  (a+3  +  8)s2  +  (8a+43+Y+16)s  +  15a  +  4y  +  4  =  0 ,   (5-14) 
with 

(5-15) 


B,    =    -14 

Cl    = 

1 

Dl   =    -4 

Fl    - 

3 

B2    =    -    7 

C2    = 

-3 

D2    =    -1 

F2    = 

-8 

(5-16) 

Substitution  into  equations  5-12  and  5-8  gives  the  singular 
line 

14a  +  133  +  35  =  0,  (5-17) 

produced  when  the  system  is  feedback  compensated  by 

Y(a,3)  -  -  |ia  -  £1*  -  |i  .  (5-18) 

The  characteristic  equation  is  then 

s3  +  (a+3+8)s2  +  (Jia  +  |p3  +  jp)  s  -  6a  -  ^3  -  §^  =  0 , 

(5-19) 

2 

which,  when  divided  by  s   +  s  +  1,  the  product  of  the  constant 

roots,  results  in  the  reduced  characteristic  equation,  or  in 
this  case,  the  remaining  root, 

s+a+3+7=0.  (5-20) 

Obviously,    for   a   stable   system, 

-a-3-7<0.  (5-21) 
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Since,  from  equation  5-17, 

14     3S 

3  =  "  rf*  "  i!  ■  (5"22) 

j-ja  -  ||  <  0  ,  (5-23) 

and 

a  <_  56  (5-24) 

for  a  stable  system  operated  on  the  singular  line. 

If  Routh's  stability  criterion  is  applied  to  equation 
5-19,  a  stability  area  in  the  parameter  plane  is  established 
by  the  boundaries  formed  by  curves 

a  +  6  +  8  =  0  (5-25) 

22a2  +  43a3  +  227a  +  21£2  +  2233  +  708  =  0     (5-26) 
12a  +  116  +  21  =  0.  (5-27) 

The  singular  line  intersects  the  boundary  of  equation  5-27  at 
a  =  56 . 

The  parameter  plane  plot  of  equation  5-19  is  shown  in 
figure  5.2,  where  the  real  root  lines  have  a  finite  error  in 
slope  because  an  unresolved  "bug"  in  the  computer  program  for 
graphing  real  root  lines. 
Example  VIII 

Type  0-4th  Order-Cascade  Compensation-Two  Parameters  in  Plant 
Denominator 

For  the  system  in  figure  5.3a,  cascade  compensation  is 
added  as  shown  in  figure  5.3b,  resulting  in  the  characteristic 
equation 
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Scale :a-50/in. ,6-50/in. 
Fig.  5.2  Example  Vll-Characteriatic  Equation 
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s5  +  (a+3+k2k1+3)s4  +  [  (k2k1+3)  (a+6)  +  ag  +  3k1k2]s3 

+  [ (3k2k1+2) (a+e)  +  (k2k1+3) gg] s2  +  [2k2k1(a+6)     (5-28) 

+  (3k  k  +2)a6+y]s  +  2k2k  aB  +  k.,Y  =  0  , 
giving,  for  roots  y   =    0 . 5 ,  oo   =1, 

Bl  =  Cl  =  2k2kl~1   Dl  =  "k2kl+2    El  =  ~kl   Fl  =  "3k2k2+1 

(5-29) 

B2  =  C2  =  k2k,+l    D   =  -2k2k1+l   E2  =  -1    Fo  =  ~k2kl"2* 

(5-30) 

For  the  case  of  the  a3  product,  equation  5-12  can  be  extended 

to 

(E2B,-E1B2)a  +  (E2C1-E1C2)e  +  (E^-j-E.^ )  a3  +  E2Fi  "  EiF2  =  °' 

(5-31) 

which  leads  to  the  singular  line 

(l-2k2k1+k2k2+k1) (a+B)  +  (k2k1-2-2k2k2+k1) a3 

+  3k2k   -  1  -  k2k2  -  2k1  =  0.  (5-32) 

For  a  linear  singular  line, 

k2kl  "  2  "  2k2kl  +  kl  =  °'  (5-33) 


which,  when  solved  for  k. ,  yields 


k„+l 


kl  "  4%"  ±   4^  /k2"14k2  +  1  •  (5"34) 


Setting  the  radical  equal  to  zero  and  solving  for  k- 


13.93  <   k2  <    0.07  (5-35) 


for  real  values  of  k, .   Choosing  k-  =  13.93  in  order  to  obtain 
the  greatest  possible  residue  for  the  singular  roots  to  in- 
crease their  dominance  makes  k,  =  0.268.   The  singular  line 
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for  the  compensated  system  is  then 

a  +  3  -  1.67  =  0.  (5-36) 

For  the  compensator,  equation  5-8  can  be  extended  to 

1   Bl    B2       Cl    C2       Dl    D2       Fl    F? 
Y(0rB,  =  -  1[(^+  Ei)a+  (_1+  _1)6  +  (gl+  gi)a6  +  gi+  ^], 

(5-37) 
which  for  k,  =  0.268  and  k?  =  13.93  gives 

_  ,  .    [14.47(a+3)-6.49a3-21.9] (s+0.268)       ,_  _0. 
Gc(s) s+3.73 (5"38) 

Although  the  pole- zero  configuration  of  the  above  comoen- 
sating  filter  is  recognized  as  undesirable  in  that  it  spans 
more  than  one  decade,  it  was  chosen  only  to  keep  the  singular 
line  linear  for  illustrative  purposes.   A  more  physically 
feasible  filter  will  result  in  a  non-linear  singular  curve. 

The  resultant  characteristic  equation  then  becomes 

s5  +  (a+3  +  6.73)s4  +  [  6 . 73 (a+3) +a3+ll . 2 ] s3  +  [13 . 2  (a+3) +6 . 73a3 ] s2 

+  [21.94(a+3)+6.71a3-21.9]s  +  3.87(a+3)  +  5.73a3  -  5.87  =  0. 

(5-39) 
The  reduced  characteristic  equation  is  then 

s3  +  (a+3+5.73)s2  +  [5 . 73 (a+3) +a3+4 . 47] s  +  6.47(a+3) 

+  5.73a3  -  10.2  =  0.  (5-40) 
Since,  from  equation  5-36, 

3  =  -a  +  1.67,  (5-41) 

equation  5-40  becomes 

s3  +  7.4s2  +  (-a2+1.67a+14.13)s  -  5.73a2  +  14.13a  +  0.6  =  0, 

(5-42) 
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for  which  a  Routh's  stability  analysis  yields  a  range  of  ooera- 
tion 

-0.04  <    a    <    2.51.  (5-43) 

Figure  5.4  shows  the  parameter  plane  plot  of  the  reduced 
characteristic  equation.   For  clarity  figures  5.5  and  5.6  show 
only  the  lines  X,   =  0.8  and  C  =  0.9,  respectively.   From  these 
plots  it  is  clear  that  the  second  pair  of  complex  roots  shift 
continuously  for  extremely  small  changes  in  a  and  3. 

Figure  5.7  is  the  parameter  plane  plot  of  characteristic 
equation  5-39. 
Example  IX 
Type  l-2nd  Order-Feedback  Compensation-No  Parameters  in  Plant 

For  the  uncompensated  system  of  figure  5.8a,  the  two 
parameters  must  be  introduced  as  in  figure  5.8b.   Tachometer 
feedback  is  used  in  order  to  retain  the  system  type.   Next 
the  singular  line  compensating  function  is  introduced  with 
accelerometer  feedback  so  that  the  singular  line  coefficients 
will  be  non-zero.   The  characteristic  equation  is  then 

(lOy+Ds2  +  (lOa+106+Ds  +  10  =  0,         (5-44) 

with,  for  L,   =    0 . 5 ,  go   =1, 

'   n 

B   =  C,  =  0     D1    =    10    F,  =  -9  (5-45) 

B2  =  C2  =  -10   D2  =  10    F2  =   0  .         (5-46) 

This  produces  a  singular  line 

10a  +  106  -  9  =  0,  (5-47) 

and 

Y(a,3)  =  *r°i  +  2-B  +  2TT  .  (5-48) 
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Fig.  5.4   Example  VIII-Reduced  Characteristic  Equation 
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Scale:a=4/in. ,3=4/in. 
Fig.    5.5      Example  VHI-Zeta  =0.8 
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Jo3 


Scale :a=4/in. ,0=4/in. 
Fig.    5.6      Example  VHI-Zeta 


=    0.9 
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Scale :a*4/in. ,0-4/in. 
Fig.  5.7  Ixample  VIII-Cb*racteristic  Equation 
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By  substitution  of  the  singular  line  relationship  into  equation 
5-48, 

Y  -  iJ  •  (5-49) 

Example  X 

Type  l-3rd  Order-Cascade  Compensation-One  Parameter  in  Plant 

Denominator 

For  the  uncompensated  system  of  figure  5.9a,  the  second 
parameter  is  introduced  via  tachometer  feedback  both  to  pre- 
serve the  system  type  and  to  avoid  a  3y  product  when  the  cas- 
cade compensator  is  added.   This  is  shown  in  figure  5.9b. 
Then  the  compensator  is  introduced  at  the  k  =  0  and  k  =  1 
levels  as  in  figure  5.9c,  yielding  the  characteristic  equation 

s   +  (a+l+k2k1)s   +  [  (l+k2k;.)  a+3+k2k  ]s   +  [k2k1  (a+$)  +4y]  s 


+  4kxy  =  0.       (5-50) 


As  usual,  for  ?  =  0 . 5 ,  co   =  1, 


Bl  =  k2kl    Cl  =  1  Dl  =  "4k2    Fl  =  _1  (5-51) 

B2  =  1      C2  =  -k2k1+l   D2  =  -4     F2  =  k2k1-l. (5-52) 

The  singular  line  is  then 

(-4k2k1+4k1)a  +  (-4-4k2k^+4kx) 3  +  4  +  4k2k^  -  Ak±  =  0.(5-53) 

If  kp  =  10  for  a  one  decade  filter  pole-zero  spread,  equation 
5-53  becomes 

-36kxa  +  (-40k^+4k1-4) 3  +  40k^  -  4k^  +4=0.  (5-54) 

The  compensating  function  is  then 

Y(a,3)  =  ^-[llk1a+(l+k1-10k^)3-l+10k^-k1]  ,  (5-55) 
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Fig.  5.9  Example  X-Type  1,  3rd  Order  Plant 
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which   gives    a    characteristic   equation 

s4    +    (a+1+lOk^s3   +    [  (l+10k1)a+$+10k1]s2 

+    [(10^  +   1J)  a  +    (5kl  +   jjL  +   1)0   -  jjL  +   5kl   -   l]s 

llkl      1    kl      2      1      2    kl 
+  -y^a  +  (j  +  J1-  ~    5kp0  -  j   +  5k^  -  ji  =  0    (5-56) 

and  a  reduced  characteristic  equation 

s2  +  (a+lOk^s  +  10k1a  +3-1=0.         (5-57) 

If  k,  =  2,  which  gives  good  dominance  characteristics  to  the 
singular  roots,  equation  5-57  becomes 

s2  +  (a+20)s  +  20a  +3-1=0;  (5-58) 
and  the  singular  line  becomes 

6a  +  133  -  13  =  0.  (5-59) 

Substitution  of  equation  5-59  into  equation  5-58  gives 

s2  +  (a+20)s  +  2-^a  =  0,  (5-60) 

which  may  be  rewritten 

,     254s 
a  ( s  +  — y^> 

s(s+20) 1  (5-61) 

for  root  locus  analysis  as  in  figure  5.10.  Obviously  the  two 
non-singular  roots  are  real,  and  the  system  is  stable  for  all 
positive  values  of  a. 

The  characteristic  equation  is 

s4  +  (a+21)s3  +  (21a+3+20)s2  +  (^-a  +  i|-B  +  ^j)  s 

+  11a  -  ije  +  2j   =  0     (5-62) 
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for  the  system  compensated  by 

(Ha  -   ZU   +   11)   (s+2) 
Vs>  "  —    's+20  16 •  <5"63> 

Parameter  plane  plots  of  the  reduced  characteristic  equation 

and  the  characteristic  equation  are  shown  in  figures  5.11  and 

5.12. 

Example  XI 

Type  l-4th  Order-Parameter  in  Numerator 

When  a  parameter  appears  in  the  numerator  of  an  uncom- 
pensated system,  the  technique  used  in  the  previous  examples 
must  be  modified  to  avoid  an  ay  or  3y  product.   In  chapter 
four  it  was  shown  that  the  plant  numerator  could  be  cancelled 
by  the  compensator,  but  in  many  situations  this  may  not  be 
desirable  or  even  possible.   However,  there  do  exist  systems, 
especially  electrical  or  electronic,  whose  transfer  function 
factors  may  be  physically  separable,  so  that  the  numerator 
parameter  may  be  removed  as  a  gain  or  filter  component  as 
shown  in  figures  5.13a  and  5.13b.   In  this  case  the  compensat- 
ing function  may  then  be  introduced  in  the  numerator  of  a  oar- 
allel  filter  having  the  same  denominator,  as  shown  in  figure 
5.13c,  where  the  compensated  characteristic  equation  is 

s4  +  (a+2)s3  +  (2a+100)s2  +  (lOOa+e+y) s  +  2g  +  0 . ly  =  0 

(5-64) 

if  k  =  0.1,  a  completely  random  value.   Then  for  £  =  0.5, 

to   =  1 , 
n 


B1  =  1      C1  =  -2    D1  =  -0.1    F1   =  98       (5-65) 
B2  =  -98    C2  =  -1    D2  =  -1      F2  =  99.      (5-66) 
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Fig.  5.11  Example  X-Reduced  Characteristic  Equation 
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Fig.  5.12  Example  X-Characteristic  Equation 
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Fig.  5.13  Example  XI -Type  1,  4th  Order  Plant 
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The  singular  line  is 

10.8a  -  1.93  +  88.1  =  0.  (5-67) 

The  compensating  function 

y(a,3)  =  -44a  -  10.53  +  539.5  (5-68) 

results  in  a  characteristic  equation 

s4  +  (a+2)s3  +  (2a+100)s2  +  ( 56a-9 . 53+539 . 5) s  -  4.4a 

+  0.953  +  53.95  =  0     (5-69) 
and  a  reduced  characteristic  equation 

s2  +  (a+l)s  +  a  +  98  =  0  (5-70) 

which  has  left-half  plane  roots  for  a  >  -1  and  real  roots  for 
a  >  20. 8. 

A  root  locus  plot  for  a  >  0  is  shown  in  figure  5.14,  and 
a  parameter  plane  plot  of  the  characteristic  equation  is  shown 
in  figure  5.15. 
Example  XII 
Type  2-6th  Order-Parameter  in  Plant  Numerator 

For  the  system  of  figure  5.16a  again  assume  it  is  pos- 
sible to  separate  the  numerator  parameter  as  in  figure  5.16b. 
Then  a  compensation  scheme  could  be  as  depicted  in  figure 
5.16c,  where  if  k,  =10,  k2  =  100,  the  characteristic  equation 

s7  +  40s6  +  (a+600)s5  +  (30a+4000)s4  +  ( 300a+3+Y+10000) s3 

+  (1000a+203+1010y)s2  +  (1503+10050y) s 

+  5003  +  50000y  =  0  (5-71) 

gives,  for  €,   =  0 . 5 ,  co   =1, 

B±    =  701    C1  =  -481    D±    =    -48991    F±   =  -9440    (5-72) 

B2  =  971    C2  =  -130    D2  =  -9040     F2  =  -3401.   (5-73) 
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Fig.  5.14   Example  XI-Reduced  Characteristic  Equation  Root  Locus 
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Fig.  5.15  Example  Xl-Characteristic  Equation 
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Fig.  5.16   Example  XII-Type  2,  6th  Order  Plant 
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The  singular  line  is 

23.2a  -  $  -  40.2  =  0.  (5-74) 

Y(a,B)  =  0.06091a  -  0.0120553  -  0.28465.      (5-75) 
The  resultant  characteristic  equation  is 

s7  +  40s6  +  (a+600)s5  +  (30a+4000)s4  +  ( 300 . 06a+0 . 98793 

+10000)s3  +  (1061. 5a+7. 82443-287. 5)s2 
+  (612.15a+28.8473-2861)s  +  3045.5a  -  102.83 
-  14230  =  0;  (5-76) 

and  the  reduced  characteristic  equation  is 

s5  +  39s4  +  (a+560)s3  +  (29a+3401)s2  +  (270 . 06a+0 . 98793 

+6039)s  +  761.44a  +  6.83653  -  9728  =  0,        (5-77) 
which  becomes,  after  application  of  the  singular  line  relation- 
ship , 
s5  +  39s4  +  (a+560)s3  +  (29a+3401) s2  +  (293 . 03a+5999 . 3) s 

+  920.73a  -  10000  =  0,  (5-78) 

for  which  a  root  locus  plot  is  shown  in  figure  5.17.   Parameter 
plane  plots  of  the  reduced  characteristic  equation  and  the 
characteristic  equation  are  shown  in  figures  5.18  and  5.19. 

Although  the  preceding  examples  are  by  no  means  the  only 
method  for  effecting  singular  line  compensation,  they  illus- 
trate the  application  of  a  mathematically  rigorous  technique 
that  is  both  physically  feasible  and  laboriously  expedient. 
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Pig.  5.17  Example  XII-Reduced  Characteristic  Equation  Root  Locus 
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Scale  :a*=100/in.  ,0-2000/in. 
Fig.  5.18  Exaapx9  x.i -Reduced  Characteristic  Equation 
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Scale :a=100/in. ,B»2000/in. 
Pig.  5.19  Example  Xll-Characteristic  Equation 
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CHAPTER  VI 
RECOMMENDATIONS  FOR  FURTHER  INVESTIGATION 

This  thesis  in  no  manner  exhausts  the  study  of  parameter 
space  and  singularity.   Although  the  compensation  techniques 
can  be  extended  to  many  types  of  control  systems,  it  would  be 
foolish  and  presumptuous  to  assume  universal  applicability. 

Solid  state  operational  amplifiers  offer  an  attractive 
means  for  synthesizing  the  compensators  discussed  in  the  ore- 
vious  chapters,  but  examination  of  the  compensator  transfer 
functions  reveal  that  truly  self-adaptive  systems  cannot  be 
realized  by  these  means  until  a  voltage  regulated  linear 
resistance  device  can  be  achieved. 

Digital  compensating  devices  should  be  investigated  for 
systems  where  expense  is  not  an  overriding  factor. 

The  possibility  exists  that  Cadena's  method  of  folding 
the  third  co-ordinate  into  the  a-8  plane  to  obtain  a  grid  of 
y   values  may  be  repeated  n-2  times  for  n  parameter  problems. 

Finally,  a  computer  program  to  analytically  solve  the 
system 

B  a  +  C,B  +  DjY  +  E,6  +  ...  +  N..V  +  P.  =  0     (6-1) 

B2a  +  C23  +  D2y  +  E26  +  ...  +  N2v  4-  P2  =  0     (6-2) 

by  iterative  processes  should  be  possible. 
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APPENDIX  I 


Table  of  Values  of  Ufc  (0.5) 


U   (0.5) 


-1  -1 

0  0 

1  1 

2  1 

3  0 

4  -1 

5  -1 

6  0 

7  1 

8  1 

9  0 
10  -1 
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